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Abstract 

We investigate channel equalization for Rayleigh fading channels under bounded 
channel uncertainties. We analyze three robust methods to estimate an un- 
known signal transmitted through a Rayleigh fading channel, where we avoid 
directly tuning the equalizer parameters to the available inaccurate channel 
information. These methods are based on minimizing certain mean-square 
error criteria that incorporate the channel uncertainties into the problem 
formulations. We present closed-form solutions to the channel equalization 
problems for each method and for both zero mean and nonzero mean sig- 
nals. We illustrate the performances of the equalization methods through 
simulations. 

Keywords: Channel equalization, Rayleigh fading, minimax, minimin, 
minimax regret. 



1. Introduction 

The problem of estimating an unknown desired signal transmitted through 
an unknown linear time-invariant channel is heavily investigated in the sig- 
nal processing literature [H, El El III El E] • Although the underlying channel 
impulse response is not known exactly in general, an estimate and an uncer- 
tainty bound on it are given [31 IH E]- In this paper we study the channel 
equalization problem for Rayleigh fading channels under bounded channel 
uncertainties [7j and investigate three robust methods to equalize Rayleigh 



* Corresponding Author. 
Email addresses: medonmez@ku.edu.tr (Mehmet A. Donmez), hinan@ku.edu.tr 
(Huseyin A. Inan), skozat@ku.edu.tr (Suleyman S. Kozat) 



Preprint submitted to Digital Signal Processing 



October 1, 2012 



fading channels. We avoid directly tuning the equalizer parameters to the 
inaccurate channel information that is available. These channel equalization 
frameworks we investigate are based on minimizing certain mean-square error 
(MSE) criteria, which incorporate the channel uncertainties into the problem 
formulations. The first approach we investigate is the affine minimax equal- 
ization method [H El E], which minimizes the estimation error for the worst 
case channel perturbation. The second approach we study is the affine min- 
imin equalization method [3 HO] , which minimizes the estimation error for 
the most favorable perturbation. The third approach is the affine minimax 
regret equalization method SI [HI E] , which minimizes a certain "regret" as 
defined in Section [2] and further detailed in Section [3} We provide closed-form 
solutions to the affine minimax equalization, the minimin equalization and 
the minimax regret equalization problems for both zero mean and nonzero 
mean signals. Note that the nonzero mean signals frequently appear in it- 
erative equalization applications [HI [12] and equalization with these signals 
under channel uncertainties is particularly important and challenging. 

When there are uncertainties in the channel coefficients, one of the preva- 
lent approaches to find a robust solution to the equalization problem is the 
minimax equalization method [91 IH E] . In this approach, affine equalizer co- 
efficients are chosen to minimize the MSE with respect to the worst possible 
channel in the uncertainty bounds. We emphasize that although the mini- 
max equalization framework has been introduced in the context of statistical 
signal processing literature [9j HI E], our analysis significantly differs since 
we provide a closed-form solution to the minimax equalization problem for 
Rayleigh fading channels. In [3j, the uncertainty is in the noise covariance 
matrix and the channel coefficients are assumed to be perfectly known. Fur- 
thermore, note that in |8], the minimax estimator is formulated as a solution 
to a semidefinite programming (SDP) problem, unlike in here. In this paper, 
the uncertainty is in the channel impulse response and we provide a explicit 
solution to the minimax channel equalization problem. 

Although the minimax equalization method is able to minimize the es- 
timation error for the worst case channel perturbation, however, it usually 
provides unsatisfactory results on the average [5]. An alternative approach to 
the channel equalization problem is the minimin equalization method [51 [10] . 
In this approach, equalizer parameters are selected to minimize the MSE 
with respect to the most favorable channel over the set of allowed pertur- 
bations. Although the minimin approach has been studied in the literature 
[HI [TOJ, however, we emphasize that to the best of our knowledge, this is the 



2 



first closed-form solution to the minimin channel equalization problem for 
Rayleigh fading channels. 

The minimin approach is highly optimistic, which could yield unsatis- 
factory results, when the difference between the underlying channel impulse 
response and the most favorable channel impulse response is relatively high 
[5]. In order to preserve robustness and counterbalance the conservative na- 
ture of the minimax approach, the minimax regret approaches have been 
introduced in the signal processing literature [31 [13J, [6] . In this approach, a 
relative performance measure, i.e., "regret", is defined as the difference be- 
tween the MSE of an affine equalizer and the MSE of the affine minimum 
MSE (MMSE) equalizer [5]. The minimax regret channel equalizer seeks an 
equalizer that minimizes this regret with respect to the worst possible chan- 
nel in the uncertainty region. Although this approach has been investigated 
before, the minimax regret estimator is formulated as a solution to an SDP 
problem [3], unlike here. In this paper, we explicitly provide the equalizer 
coefficients and the estimate of the desired signal. 

Our main contributions are as follows. We first formulate the affine equal- 
ization problem for Rayleigh fading channels under bounded channel uncer- 
tainties. We then investigate three robust approaches; affine minimax equal- 
ization, affine minimin equalization, and affine minimax regret equalization 
for both zero mean and non-zero mean signals. The equalizer coefficients, 
and hence, the MSE of each methods have been explicitly provided, unlike 
in [HI HI El El M- 

The paper is organized as follows. In Section |2j the basic transmission 
system is described, along with the notation used in this paper. We present 
the affine equalization approaches in Section [3] First, we study the affine 
minimax equalization tuned to the worst possible channel filter. We then 
investigate the minimin approach and the minimax regret approach, and 
provide the explicit solutions to the corresponding optimization problems. 
In addition, we present and compare the MSE performances of all robust 
affine equalization methods in Section |4j Finally, we conclude the paper 
with certain remarks in Section [5] 

2. System Description 

In this section, we provide the basic description of the system studied in 
this paper. Here, the signal x t is transmitted through a discrete-time time- 
varying channel with an impulse response h t . The received signal y t is given 
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by 

y t = x t h t + n t , (1) 

where the observation noise n t is independent and identically distributed 
(i.i.d.) with zero mean and variance a\ and independent from xt- However, 
instead of the true channel impulse response, an estimate of h t is provided 

as ht, where Sh t = h t — h t is the uncertainty in the channel impulse response 
and is modeled by \h t — h t \ = \5h t \ < e, e > 0, e < oo, where e or a bound on 
e is known. 

We then use the received signal y t to estimate the transmitted signal x t 
as shown in Fig. 1. The estimate of the desired signal is given by 



x t = w t y t + k 

= w t (x t h t + n t ) + l t , 



(2) 



where w t is the equalizer impulse response. We note that in ([2]), the equalizer 
is "affine" where there is a bias term l t since the transmitted signal Xt, and 
consequently the received signal yt, are not necessarily zero mean and the 

mean sequence yt = E[y t ] is not known due to uncertainty in the channel. 

Even under the channel uncertainties, the equalizer impulse response Wt 
and the bias term l t can be simply optimized to minimize the MSE for the 
channel that is tuned to the estimate h t , which is also known as the MMSE 
estimator [T4]. The corresponding impulse response and the bias term are 
given by 



{wo,t,k,t} = argminE 

w,l 



X t 



w 



h t x t + Tit) - I 



However, the estimate 
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x 0) t = Wo,tVt + lo,t 



may not perform well when the error in the estimate of the channel impulse 
response is relatively high [l"5t |3], Hj. One alternative approach to find a 
robust solution to this problem is to minimize a worst case MSE, which is 
known as the minimax criterion, as 



{wi,uh,t} — argmin max E 

w,l \8h t \<e 



Xt-w((ht + 5h t )x t + n t ) - I 
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where wxj and minimize the worst case error in the uncertainty region. 
However, this approach may yield highly conservative results, since the esti- 
mate 



xi,t = wxtyt + I 



is formed by using the equalizer impulse response Wi jt and the bias term li t 
that minimize the worst case error, i.e., the error under the worst possible 
channel impulse response |3J Hj. Instead of this conservative approach, 
another useful method to estimate the desired signal is the minimin approach, 
where the equalizer impulse response and the bias term are given by 



Xt — w((ht + 5h t )x t + n t ) - I 



{ w 2,t, h,t} — argmin min E 

w,l \5ht\<e 

where W2,t and l-if, minimize the MSE in the most favorable case, i.e., the 
MSE under the best possible channel impulse response [5]. The estimate of 
the transmitted signal x% is given by 

A , 

x 2 ,t = W2,tVt + h,t- 

A major drawback of the minimin approach is that it is a highly optimistic 
technique, which could yield unsatisfactory results, when the difference be- 
tween the actual and the best channel impulse responses is relatively high 
0- 

In order to reduce the conservative characteristic of the minimax approach 
as well as to maintain robustness, the minimax regret approach is introduced, 
which provides a trade-off between performance and robustness [31 [HI [B] . In 
this approach, the equalizer impulse response and the bias term are chosen 
in order to minimize the worst-case "regret" , where the regret for not using 
the MMSE is defined as the difference between the MSE of the estimator and 
the MSE of the MMSE, i.e., 

{w 3 ,t,h,3} = argmin max 

w,l \8ht\<e 



E 



x t -w((h t + 5h t )x t +n t ) - I 



mini? 

w.l 



x t - w \ {h t + Sh t )x t + n t ) - I 



The corresponding estimate of the desired signal x% is given by 



A 



%3,t = w 3)t y t + l 3>t . 
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In the next section, we investigate and provide closed form solutions for 
the three equalization formulations: 

• Affine minimax equalization framework, 

• Affine minimin equalization framework, 

• Affine minimax regret equalization framework. 

We first solve the corresponding optimization problems and obtain the es- 
timates of the desired signal. We next compare their mean-square error 
performances in Section |4} 



3. Equalization Frameworks 

3.1. Affine MMSE Equalization 

In this section, we present the affine MMSE equalization framework for 
completeness jTU |TJ] . Since the channel impulse response h t is not accurately 
known but estimated by h t , a linear equalizer that is matched to the estimate 
ht and minimizes the MSE can be used to estimate the transmitted signal 
x t . The corresponding equalizer impulse response Wo s t and the bias term l 0jt 
are given by 



{wo,t,k,t} = argminE 

w,l 



x t - w ( h t x t + n t ) - I 



(3) 



Since the cost function in (|3| is convex in w and I, the minimizers Wqj and 
lot are given by 



~h t ol 



w ,t 



hWx + <rl 



3.2. Affine Equalization Using a Minimax Framework 

In this section, we investigate a robust estimation framework based on 
a minimax criteria [1^ [16], HD]- We find the equalizer impulse response wij 
and the bias term li t that solve the following optimization problem: 



{ w i,t,h,t} — argmin max E 

w,l \Sht\<e 



Xt-w[(h t + 5h t )x t + n t 



(4) 
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In (HI), we seek to find an equalizer impulse response wi >t and a bias term 
lij that perform best in the worst possible scenario. This framework can be 
perceived as a two-player game problem, where one player tries to pick Wxj 
and li jt pair that minimize the MSE for a given channel uncertainty while the 
opponent pick Sh t to maximize MSE for this pair. In this sense, this problem 
is constrained since there is a limit on how large the channel uncertainty Sh t 
can be, i.e., \5h t \ < e where e or a bound on e is known. 

In the following theorem we present a closed form solution to the opti- 
mization problem Q. 

Theorem 1: Let x t , yt and n t represent the transmitted, received and noise 
signals such that y t = h t x t + n t , where h t is the unknown channel impulse 
response and n t is i.i.d. zero mean with variance a 2 . At each time t, given 



an estimate ht of ht satisfying \h t — h t \ 
problem 



< e, the solution to the optimization 



{wi,tih,t} — argmin max E 

w,l \Sht\<e 



x t - w ((h t + Sht)x t + n t ) - I 



(5) 



is given by 



[ht-e 
J_ 
ht 



: ~h t ta 2 x < t 2 a 2 x + a 2 n 
: ~h t ea 2 x > e 2 a 2 x + a 2 



and 



■>'l a n 



: h t ea 2 . < e 2 al + a 2 



: h t ea 2 x > e 2 a 2 x + a 2 n . 



Proof: Here, we find the equalizer impulse response Wij and the bias term 
li t t that solve the optimization problem in ([5]) . To accomplish this, we first 
solve the inner maximization problem and find the maximizer channel uncer- 
tainty 5h* t . We then substitute 5h* t in ^ and solve the outer minimization 
problem to find Wx t t and li >t . 

We solve the inner maximization problem as follows. We observe that 
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the cost function in Q can be written as 



E 



x t - w ( (h t + 5h t )x t + n t ) - I 



(6) 



E [x 2 t ] + w 2 (h t + ShA + I 2 - 21E [x t ] - 2w (h t + 5hA E [x 2 ] + 2wl (h t + ShA E 
w 2 (~h t + 5h t ) + 2w (~h t + 5h t ) (lx~ t - x 2 ) + G u 



where x t — E [x t ], x\ — E [x 2 ] and Ci = E [x 2 ] + w 2 a 2 n + I 2 — 2lx t does not 



depend on 5h t . If we define a = x\ > 0, b = lx t — x 2 , u = w (h t + 5h t ) and 
C 2 = Ci — then ^ can be written as 

2 



E 



x t -w ( (h t + 8ht)x t + n t ) - I 



a ( u + - ) + C 2 , 

a 



where C 2 is independent of 5h t . Hence the inner maximization problem in 
^ can be written as 



Sht = arer max E 

\8h t \<e 



x t -w[(h t + 5h t )x t + n t ) - I 



arg max \w\ 

\Sht\<e 



Sh t + 
(7) 



lx t - 

x t 



vox; 



If we apply the triangular inequality to the second term in (|7|), then we get 
the following upper bound: 



\w\ 



Sfh + h t + 



wxt 



— \ w \ 


\Sth\ + 


h t + — 




— \ w \ 


e + 


~h t + lXt 








1 

WXt 








wx^ 





where the upper bound is achieved at Sh* = esgn ( h t + xt J^* \ ( where sgn(z) 

V wx t 

1 if z > and sgn(z) = — 1 if z < 0. Hence it follows that 



8h*+ = arg max E 

\Sht\<e 



x t -w[(h t + Sh t )x t + n t ) - I 



h t + lxt - X < > 0, 



wxt 



/ <~v> nn 2 

-e :h t + — 1 < 0. 
wx 2 

(8) 
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We next solve the outer minimization problem as follows. We first note 
that the minimum in ^ is taken over all w G R and I G R. If we write u = 

[w,l] G R 2 in a vector form and define U = < u = [w, iy g m 2 1 h t + 



tot x t ^> Q 



U 



\w, if G R 2 



ht + 



lxt—x? 



< 



\U, i.e. 



then it follows that V 

W U V = M 2 and U fl V = 0. Hence, the cost function in the outer minimiza- 
tion problem in ^ is given by 



max E 

\5h t \<e 




x t -w((h t + 6h t )x t + n t ) -I 



x t — w \{h t + e)x t + n t ) — / 



x t - w ( (h t - e)x t + n t ) -I 



: [w,lf EU, 
: [w, if G V 



We first substitute Sh t = e and find the corresponding {w, 1} pair that min- 
imizes the objective function in ^ to check whether [w, I] G U. We next 
substitute 5h t = — e and find the corresponding {w, 1} to check whether 
[w, I] G V. Based on these criteria, we obtain the corresponding equalizer 
impulse response and the bias term pair {wij, l\,t\- 

We first substitute Sh t = e in the objective function of ^ to get the 
following minimization problem: 



{w*, I*} = arg min 



w.l 



(9) 



x 2 t + w 2 [ ( h t + e ) xf + a 2 n ) + I - 2/^ + 2tuZ ( ^ + e ) ^ - 2w ( h t + ej x\ 

(10) 



We observe that the cost function in (10) is a convex function of w and I 
yielding 



w ' 



h t + e)a 



h t + e) al + al 



cri 



where a 2 — E Ux t — x t ) 2 ] . However we have 



xr 



l*xi 



■urx 



{ht + e)al 



> h. 



(11) 
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so that [w*,l*f <£U. 

We next substitute 5h t = — e in the cost function of ^ to get 

{w*,l*} = argmin (12) 

w,l 

jx? + w 2 (jfi t - x\ + o-^J + / 2 - + 2wZ (fi t -e^x~ t -2w (h t - x\ 

(13) 



The cost function in (13) is also a convex function of w and I so that we get 



W = 7^ ^ , I = 

ht-e) a 2 x + a 2 n (h t - e) o 2 x + o 2 n 



If the condition h t ea x < e a x + a„ holds, then we have 

7 7 °n X + ~~ l Xi 

h t <h t -e + -z =-= < 



{ht — e)x 2 wx\ 

so that [w*,l*] T G V. Thus, the corresponding equalizer impulse response 
and the bias term are given by = (j^O^g — anc j ^ _ — — _^ 

(rat— ej o-^+it^ (rat— ej o x +(T n 

respectively. However, if the condition h t ea 2 < e 2 a 2 + a 2 does not hold, then 
[w*,l*] T ^ V so that the corresponding equalizer impulse response and the 

2 

bias term are given by are W\^ t = t~ and l\ t = ,r — 2- Hence, we have 

Wi,t = \ [ h t- e ) , 
l t : ~h t ea 2 x > t 2 a 2 x + a 2 n 

■ ~ h t^l < t 2 °l + °l 
: h t ea 2 > e 2 a 2 + a 2 n 



l lt = { {h-efvl+vl 



Hence, the proof follows. ■ 

In the following corollary, we provide a special case of Theorem 1, where 
the desired signal xt is zero mean. 

Corollary 1: When the transmitted signal xt is zero mean, the solution to 



10 



the optimization problem 



{ w i,tih,t} — argmin max E 

w,l \6ht\<e 

is given by 



(x t - W ((h t + Sh t )x t + n t )-l 



(14) 



ht-e) 



ht 



,a' i ■.e[h t -e)<\ 



r i ht - e I > | 



A 



where 5* = cr^/a^ is the signal-to-noise ratio (SNR). 
Proof: The proof directly follows from Theorem 1, therefore, is omitted. 



3.3. Affine Equalization Using a Minimin Framework 

In this section, we study the minimin equalization framework, where the 
inner maximization of the minimax framework is replaced with a minimiza- 
tion over the uncertainty set [5J [T71 [TO]. We seek to solve the following 
optimization problem: 



{w2,tih,t} = argmin min E 

w,l \8ht\<e 



x t -w[(h t + Sh t )x t + n t ) - I 



(15) 



The following lemma is introduced to demonstrate that min operators in 



(15) can be interchanged, which will be used in Theorem 2. 



Lemma 1: For an arbitrary function f(x, y, z) and nonempty sets X, y and 
Z, we have 

min min f (x, y, z) = min min f(x,y,z), 

assuming that all minimums are achieved on the corresponding sets. 
Proof: The proof is straightforward hence omitted here. 

In the following theorem we present a closed form solution to the opti- 



mization problem ( 15 ) 



Theorem 2: Let Xt, yt and nt represent the transmitted, received and noise 
signals such that y t = h t x t + n t , where h t is the unknown channel impulse 
response and n t is i.i.d. zero mean with variance o\. At each time t, given 
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an estimate h t of h t satisfying \h t — h t \ < e, the solution to the optimization 
problem 



{w2 ; t,h,t} = argmin min E 

w,l \Sht\<e 

is given by 



w 2 ,t 



x t - w I (h t + 5h t )x t + n t ) - / 



h t + esign(^) ) a 2 x 



(16) 



and 



h t + esign(^) ) a 2 x + a\ 



x t crt 



2,t 



ht + esign(/i t )J a2 + a\ 
Proof: Here, we find the equalizer impulse response w 2 1 and the bias term 



l 2 t that solve the optimization problem in ( 18 ). We first note that, by Lemma 



1, we can interchange min operators in (18 
in (18) is equivalent to 



so that the optimization problem 



min min E 

w,l \Sht\<e 



x t -w((h t + 5h t )x t + n t 



min mini? 

\5ht\<e w,l 



x t — w[(ht + Sh t )x t + n t ) - I 



(17) 



Hence, we first solve the inner minimization problem in (17) and find the 
minimizers w* and I*. We then substitute w* and I* in (17) and solve the outer 
minimization problem to find the minimizer 5h* t , which yields the desired 
equalizer impulse response w 2 t and l 2t - 



We observe that the objective function in (fT5|) can be written as 
E (x t -w ((ht + 5h t )x t + n^j - l) 

= E [x\] + w 2 (h t + 5h t y + l 2 - HE [x t \ - 2w (h t + 5h^j E [x 2 t ] 
+ 2wl (h t + 5h t ) E [x t ] 



9 2 

Xt + w 



(18) 



h t + Sh t x 2 t + < + V - 2lxl + 2wl [h t + 8ht)xt 



-2w[h t + 5h t ) x 2 , 
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We first solve the inner minimization problem in the right hand side of 
(17) with respect to w and I as follows. If we define F(w,l,Sh t ) 



E 



x t -w((h t + 5h t )x t + n t ) - I 



, then we set the first derivatives 



of F(w,l,8ht) with respect to w and I which yields the minimizers w* and 



OF 



and ^ 



0. The 



i=i* 



I*, respectively, i.e., w* and I* satisfy dw 

corresponding partial derivative of the cost function F(w, I, 5h t ) with respect 
to I is given by 

dF 



i=i* 



21* - 2x1 + 2w* [ht + Sh t )xl = 



dl 

so that I* = x~t — w* {lit + 5htj x~t. The corresponding partial derivative of 
F(w, I, 5h t ) with respect to w is given by 



dF 

dw 



w=w* 



2w* Iht + Sht) x\ + al )+ 21* (h t + 5h t ^ - 2 [h t + 8h t )x 



which implies that w* 



(ht+Sht) x'l+al 

~h t + 5h t ) a 



. Thus, we get that 



w 



h t + 5h t ) al + at 



I* 



x t <J 2 n 



for a given Sh t . 

We next solve the outer minimization problem. If we substitute w* and 
I* in F(w, I, 5h t ), then we obtain 



5h* = arg min min E 

\8ht\<e w,l 



x t -w((h t + Sh t )x t + n t ) - I 



arg min F(w*,l*,5h t ) 

\5h t \<e 



arg mm 

|<5fct|<e 



arg max 

\6h t \<e 



ht + Sht) al + al 



h t + Sht 



(19) 
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so that 8h* t = esign(/i t ). Hence, the equalizer impulse response w 2>t and the 
bias term Zg,* are given by 



^ + esign(^) ) cr 



/i t + esign(/i t )J cr| + a\ 

2 



X t°r, 



h t + esign(/i t ) J o| + a 2 



Hence, the proof follows. ■ 

In the following corollary, we provide a special case of Theorem 1, where 
the desired signal Xt is zero mean. 

Corollary 2: When the transmitted signal xt is zero mean, the solution to 
the optimization problem 



{w2,t,h,t} — argmin min E 

w,l \Sht\<e 



x t - w ( (ht + 5h t )x t + n t ) - I 



is given by 



U>2,t 



h t + esign(/i t 



ht + esign(/i t 



+ 



and 



A 



where S = a 2 x /a 2 n is the SNR. 
Proof: The proof follows from Theorem 2 when x~l — 0. ■ 

3.4- Affine Equalization Using a Minimax Regret Framework 

In this section, we investigate the minimax regret equalization framework, 
where the performance of an affine equalizer is defined with respect to the 
MMSE affine equalizer that is tuned to the unknown channel [31 El [Til EH] ■ 
We emphasize that the minimax equalization framework investigated in Sec- 



tion 3.2 may produce highly conservative results since the equalizer impulse 



response w and the bias term / are optimized to minimize the worst case MSE 
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j. Moreover, the minimin equalization framework introduced in Section 3.3 
is a highly optimistic method where the equalizer parameters are optimized 
to minimize the MSE that corresponds to the most favorable channel [5]. 
Thus, the minimin approach may also yield unsatisfactory results in certain 
applications, where the channel estimate is highly erroneous [3j. In this con- 
text, the minimax regret equalization framework can be used to improve the 
equalization performance while preserving the robustness [6] . In this ap- 
proach, we find the equalizer impulse response wsj and the bias term l 3jt that 
solve the following optimization problem: 



{w 3 ,t,h,t} = argmin max 

w,l \Sht\<e 



E 



x t -w( (h t + 5h t )x t + n t ) - I 









— mini? 




w,l 



x t -w[(h t + 5h t )x t + n t 
(20) 



We note that from Section 3^, the solution to the minimization problem in 
the objective function is given by 



min E 

w,l 



x t -w[(h t + 5h t )x t + n t ) - I 



h t + 5h t ) (Jl + ol 



Hence the optimization problem in (20) is equivalent to 



argmin max 

w,l \8ht\<e 



E 



x t -w[(h t + Sh t )x t + n t ) - I 







)-')[ 


— min E 




w,l 



x t -w[(h t + 8h t )x t + n t 



argmin max \ E 

w,l \5ht\<< 



x t - w [(h t + 5h t )x t + n t ) - I 



ht + Sht) +o 2 n 
(21) 



We first expand the term 



h t + Sh t ) al + a 2 n 
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in (21) around Sh t = yielding 



h + Sh t ) + a 



hi + al 



h 2 t (?l + el 



Hence, instead of (20), we solve the following optimization problem: 



{w 3 ,t,h,t} = argmin max 

w,l \5h t \<e 



E 



x t - w \ {h t + 5h t )x t + n t ) - I 



+ Sht 



hi + ol 



hhl + °l 



(22) 



which provides satisfactory results even under large derivations 6h t as shown 
in the Simulations section. 

In the following theorem we present a closed form solution to the opti- 



mization problem (22) 



Theorem 3: Let Xt, yt and nt represent the transmitted, received and noise 
signals such that yt = htXt + n t , where ht is the unknown channel impulse 
response and n t is i.i.d. zero mean with variance a 2 . At each time t, given 
an estimate h t of h t satisfying \h t — h t \ < e, the solution to the optimization 
problem 

{w 3 ,t,h,t} = argmin max 

w,l \5h t \<€ 



E 



x t - w [{h t + Sh t )x t + n t ) - I 



h 2 t +< 



+ Sh t 



hWx + a\ 



(23) 



is given by 



3,fJ 



H, it] 
KG] 



f>0,g>0, 
f < 0, g < 0, 
/ > 0, g < 0, 
/ < 0, g > 0, 
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where 



H, it] 



h t + e I a 



x t a 2 



ht + e) a 2 + a 2 ( h t + e) o 2 x + a n 



ht-e) a 



ht-e) (T 2 x + a 2 n [ht-e] cr 2 x + a 2 n 



[u>3, 11} = arg min 

Me-Q^*],^*]} 



max < E 

\Sht\<e 



x t — w[(ht + 5h t )x t +n t ) - I 



+§ht 2h t ala: 



2^4 

x 



h 2 t +a 2 n 



[w*4, 11} = arg min < E 

[w,l] 



x t -w[ h t x t + n t ) - / 



°Wx 

M + < 



A 

9 = e 



a , 



+ 



h t a, 



((u 



h t + e)\ 2 x + a 2 (~ht + e)a 2 (ht + e" 



2 \ 

rr 

x 1 n 



h t + e) al + a 2 ^ 



ht-e) o* + o* [ht-e)a 2 x (~ ht 



+ 



ht°l 



hW x + a\ 

\ 

( (~h t -e\ o 2 x + o 2 y 



\ 



h\°l + °l 



Proof: We first observe that the objective function in (23) can be written 



17 



as 



E 



x t - w ( ( h t + 5h t ) x t + n t ) - I 



t— X 2 



h 2 + ol 

x\ + w 2 ( ht + Sh t J x 2 t + w 2 a 2 n + I 2 - 2lxl -2w [h t + Sh t ) x'f + 2wl [h t + 5h t )x~ t 
2h t a 2 n a A x 



2 2 

+ oht 



h t + < (h 2 a 2 x + a, 



( 



w 2 (ht + Sh t x\ + (ht + 5h 



2wlx t — 2wx\ + 



2VX 



V 



(24) 



A 



where Di = + wV„ + r — 2/^ 



ht , a " 7 °'' I ,i is independent of <5/i t 



If we define a = w 2 x 2 t >0,b = 2wlx- t - 2wx 2 + /r f ta } u \ 2 and D 2 = D : b * 



then (24) can be written as 

E ( x t - w ( (h t + 5h t ) x t + n t ) - I 



2ht<j\ 



4a' 



h 2 t + a 2 n 



a \ u + 



2a 



D 



2- 



where D 2 is independent of 5h t . Hence, the inner maximization problem in 



(23) is given by 



8hl 



arg max < E 

\5ht\<e * 



x t -w((ht + Sh t )x t + n t ) - I 



h 2 t + o 2 n 



arg max 

|<5ht|<£ 



Sht + h + 



lx t 



wx: 



w 



+ 



hto~p% 



w 2 x 2 h 2 a 2 x + a 2 n 



(25) 
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By applying the triangular inequality to the cost function in (25), we get the 
following upper bound: 



r Ixt 1 
Sh t + h t + ^-- + 

wx\ w 



ht<J 2 n <J 



2^4 

x 



w 2 x\ ( h 2 t a 2 + a 2 



< \Sht\ + 



< e + 



h t + 



lx t 



WX7 



lxt 



wx; 



1 

w 



+ 



w 2 x 2 { ~h 2 t a 2 x + a 2 



htala 



2^4 

x 



W 



w 2 x 2 t (h 2 t al + al 



where the upper bound is achieved at 5h% = esgn h t + — - -| — ^ 



Hence it follows that 



8K 



arg max < E 

\Sht\<e 1 



' 1 ^ Xt 

e : h t -\ = 



wx 



-e :h t + 



lx t 



wxt 



x t -w((h t + 8h t )x t + n t ) - I 



°n a l +Sht 2h t a l°x 



h 2 + a 2 



1 



w 



w 2 x\ ( h 2 a 2 + cr 2 



+ 

w 



2^4 

x 



w 2 x 2 t (h 2 al + al 



> 0, 



< 0. 



(26) 



We next solve the outer minimization problem as follows. If we write u 



then it follows that 



[w, if G M 2 and define M = {u = [w, if e M 2 \ h t + % - - + — = 



Af 



u = [w, if G 



h t 



lx t 



M U Af 



< 



\M, i.e. 



> 



I 2 and A4 PI M = 0. Hence, the cost function in the outer 
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minimization problem in (23) is given by 



max < E 

\Sh t \<e 



E 

E 



x t -w((h t + 6h t )x t + n t ) - I 



°n a l +Sht 2 ^ CT n^ 



h\ + at 



hWx + °l 



x t -w[(h t + e)x t + n t ) - I 



x t -w[(h t - e)x t + n t ) - I 



_ _|_ c g^tqjoj 



— e 



2h t crlcrt 



: [w,iy e M 

: [w, l] T G Af. 



We first substitute Sh t = e and find the corresponding {w, 1} pair that min- 



imizes the objective function in (23) to check whether [w,l] gM. We next 
substitute 5h t = — e and find the corresponding {w, 1} to check whether 
[w, 1} G M . Based on these criteria, we obtain the corresponding equalizer 
impulse response and the bias term pair {w 3jt , / 3>i }. 



We first substitute Sh t = e in the cost function in (23) to get the following 
minimization problem: 

{wl, l{} = arg min 

w,l 



x\ + w 2 ( h t + e x\ + w v n + r -2w ih t + e) x\ - 2x~ t l - 2x~ t w [h t + e)l 



"7, "7- +e _ -l>l(T„CT, 



(27) 



Since the cost function in (28) is a convex function of w and I, we get that 

(ht + e) ol 
w* = . V J , 11 = 



x t cr; : 



l - N 2 

h t + e) al + al 



ht + e) a 2 x + al 



We observe that [wl, l\] G M if and only if 



Xt 2 (T 2 n 



- + 

~h t + e\al + al [fh + e]^ (h t + e 



h\ol + ol 



> 0. 
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We next substitute 5h t = — e in the cost function in (23) to get the 
following minimization problem: 



= argmin 

w.l 



x'f + w 2 ( ht - e ) x'f + w 2 a 2 n + I 2 - 2w ( h t - e ) x'f - 2x t l - 2x t w [h t — e)l 



2VM 



h 2 + 



t T(T< 



h 2 a 2 x + a) 



(28) 



Since the cost function in (28) is a convex function of w and I, we get that 



?/>.; 



ht — e) a 



h t — e) al + at 



- 1* 



ht — e) ol + ol 



Note that [W2V2] e N if an d only if 



A 

9 = e 



(h t -e)\l + ol {ht-e)a 2 (h t 



h t at 



IK 



h t -e) al + a 2 ^ 



hW x + a* 



< 0. 



There are four cases depending on the values of h t , e, x t , x 2 , a 2 : 

• Case 1: / > and g > 0. 
In this case, we have 



w 3 ,t 



h t + e)a. 



h t + e) ol + ol 



and 



h t + e) ol + at 



since [wl,ll] G M and [wl,l* 2 ] £Af. 
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Case 2: / < and g < 0. 
In this case, we have 



U>3,t = 



(h t - e j cr 



(ht-e) al + ai 



and 



^3,t — 



x t a 2 n 



(h t -e) al + al 

since [k;*, /*] and [tug, ^] G A/". 
Case 3: / > and g < 0. 

In this case, we have ZJ] G and [u^, ^] G A/" so that 



3,tJ 



arg mm 

[w,i]e{[w*,q],[w* 2 ,i* 2 ]} 



max < 

\Sht\<e 



x t -w {[h t + 5/i t )x t + n^j - I s ) 



2 2 

h 2 t + al 



+ Sh t 



2hi 



Case 4: / < and g > 0. 

In the last case, we have the optimum points on the curve h t + -= — 

- H — h :t cr "' 7x — T2 = 0. Therefore 5h* t = and the corresponding coef- 

ficients are given as the solution to the following optimization problem: 



[w 3 ,t,h,t] = argmin <^ E 

[w,l] 

subject to 



h t + 



lXt 



+ 



(x t - w (h t x t + n^j - l) 
h t a 2 n a A x 



h 2 t + < 



wxi 



w 



= 0. 



2^2 



wx 



Hence, the proof follows. 
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4. Simulations 

We provide numerical examples in different scenarios in order to illus- 
trate the performances of the equalization methods. We first illustrate the 
performances of the channel equalization methods for a given perturbation 
bound. We demonstrate that the the minimax equalization method yields 
the best worst case MSE performance among all methods for these simula- 
tions since it optimizes the worst case MSE with respect to the worst case 
channel impulse response. We next present the average MSE performance of 
each method over different channel perturbations. We show that the mini- 
max regret method has better average MSE performance than the minimax 
and minimin equalization methods for these simulations. 

In the first set of experiments, we randomly generate a transmitted signal 
x t with mean 0.01 and variance 1. We also generate a Gaussian channel noise 
n t with zero mean and unity variance. The channel estimates are constructed 
using ht = ht + Sht, where ht = 1.05 and the perturbation 5h t is randomly 
generated from a zero mean and e standard deviation Gaussian distribution 
and truncated to give \5h t \ < e with e = 0.03 for each trial. Here, we la- 
bel the method in Theorem 1 as "Minimax", the method in Theorem 2 as 
"Minimin" , and finally the method in Theorem 3 as "Minimax regret" . For 
each method and for each random perturbation, we find the corresponding 
equalizer parameters Wt and U to calculate the estimates of the transmit- 
ted signal Xt- After we calculate the mean-square errors for each method 
and for all random perturbations, we plot the corresponding sorted errors 
in ascending order in Fig. 2 for 200 trials. Since the minimax equalization 
method optimizes the worst case MSE with respect to worst possible pertur- 
bation, it yields the smallest worst case MSE among all methods for these 
simulations. However, the overall performance of the minimax method is sig- 
nificantly inferior to the minimin and the minimax regret methods due to its 
highly conservative nature. Furthermore, we notice that the minimax regret 
method provides better average performance compared to the minimax and 
the minimin methods and superior worst case performance compared to the 
minimin method for these simulations. 

For the second experiment, we randomly generate 200 random perturba- 
tions Sh t , where \Sh t \ < e for different perturbation bounds and compute the 
averaged MSEs over 200 trials for the minimax, minimin and the minimax 
regret methods. In this case, we randomly generate a transmitted signal Xt 
with zero and variance 1. The channel noise rit is generated from a Gaussian 
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distribution with zero mean and unity variance. Here, we construct the esti- 
mates of the channel impulse response by h t = ht + Sht, where ht = 1.05 and 
the perturbation Sht is randomly generated from a zero mean and e standard 
deviation Gaussian distribution and truncated to give \5h t \ < e. In Fig. 3, we 
present the averaged MSEs for each method where the perturbation bound 
varies, e G [0.1,0.3]. We observe that the minimax regret method has the 
best average MSE performance over different perturbation bounds compared 
to the minimax and the minimin equalization methods. 

In this section, we presented the performances of the minimax, minimin 
and minimax regret channel equalization methods through simulations. We 
demonstrated that the minimax approach leads to a better worst case MSE 
performances than the minimin and minimax regret approaches for these 
simulations. We also presented the average MSE performances of the equal- 
ization methods over different channel perturbations and showed that the 
minimax regret equalization method has the best average MSE performance 
among all methods for these simulations. 

5. Conclusion 

In this paper, we investigated the channel equalization problem for Rayleigh 
fading channels when the channel impulse response is not accurately known. 
We analyzed three robust methods to equalize Rayleigh fading channels that 
incorporate the channel uncertainties into the problem formulation. We first 
studied the afline minimax channel equalization framework that optimizes 
equalizer parameters to minimize the worst case MSE in the uncertainty re- 
gion. We next investigated the afline minimin channel equalization method, 
which minimizes the MSE for the most favorable channel impulse response 
in the perturbation bounds. Finally, we analyzed the afline minimax regret 
channel equalization framework, which minimizes the worst case regret in 
the uncertainty region. We explicitly provide the equalizer coefficients and 
the estimates of the desired signal for each method and for both zero mean 
and nonzero signals. We illustrated the performances of these equalization 
methods through simulations. 
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Figure Captions 



Figure 1. A Basic Affine Equalizer Framework 

Figure 2. Sorted MSEs for the minimax, minimin and minimax regret equal- 
ization methods over 200 trials when e = 0.3. 

Figure 3. Averaged MSEs for the minimax, minimin and minimax regret 
equalization methods over 200 trials when e G [0.1,0.3] 
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Figures 
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Figure 1: A Basic Affine Equalizer Framework 
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Figure 2: Sorted MSEs for the minimax, minimin and minimax regret equal- 
ization methods over 200 trials when e = 0.3. 



28 




20 40 60 80 100 120 140 160 180 200 

trials 



Figure 3: Averaged MSEs for the minimax, minimin and minimax regret 
equalization methods over 200 trials when e G [0.1, 0.3] 
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